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Measuring the hyperfine structure (HFS) of long-lived 3 P2 states of divalent atoms may offer the 
opportunity of extracting relatively unexplored nuclear magnetic octupole and electric hexadecapole 
moments. Here, using relativistic many-body methods of atomic structure and the nuclear shell 
model, we evaluate the effect of these higher nuclear moments on the hyperfine structure. We find 
that the sensitivity of HFS interval measurements in 87 Sr needed to reveal the perturbation caused 
by the nuclear octupole moment is on the order of kHz. Results of similar analyses for 9 Be, 25 Mg, 
and 43 Ca are also reported. 
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I. INTRODUCTION 

Progress in cooling and trapping techniques has en- 
abled a number of advances in precision measurements of 
atomic properties. Compared to more traditional beam 
spectroscopy, trapping atoms or ions increases interro- 
gation times and enhances spectral resolution. This is 
a well-recognized experimental trend spanning from the 
searches for permanent electric-dipole moments to atomic 
clocks. Here, we point out that measuring the hyperfine 
structure of long-lived 3 P2 states of cold divalent atoms 
may facilitate extracting the so far relatively unexplored 
nuclear magnetic octupole moments. 

A nucleus may be described as a collection of elec- 
tromagnetic moments. Generally, a nucleus of spin 
I possesses 21 moments beyond the monopole. Most 
well-known are the magnetic dipole (/x) and electric 
quadrupole (Q) moments. These two leading moments 
have been studied extensively and have contributed to 
the understanding of nuclear forces. For example, dis- 
covery of the quadrupole moment of the deuteron led to 
introducing nuclear tensor forces Here, we are inter- 
ested in the next two moments in the multipolar hierar- 
chy: nuclear magnetic octupole (Q) and electric hexade- 
capole (n) moments. 

Electromagnetic fields of various nuclear moments per- 
turb the motion of atomic electrons and reveal themselves 
in the hyperfine structure (HFS) of atomic and molec- 
ular levels. The resulting hyperfine intervals are con- 
ventionally parameterized in terms of the HFS constants 
A,B,C,D,- ■ ■ , each being proportional to the relevant 
nuclear multipolar moment. The perturbation becomes 
weaker as the rank of the multipole moment increases, 
and detecting the influence of higher multipole moments 
on the HFS intervals requires increasingly longer interro- 
gation times. Also, to observe the effects of higher rank 
multipoles on the HFS structure one is forced to work 
with electronic states of relatively high angular momenta 



J since the 2 fe -pole moment only contributes in first order 
to the HFS structure of levels with J > k/2. These two 
factors have limited the determination of the octupole 
HFS constants to a small number of isotopes of CI 0], 
Ga 3, Br In @], V @, Eu Q, Lu @, and Hf @. All 
of the above measurements have been carried out cither 
on the ground state or on metastable states. The most 
recent octupole moment measurement was carried on the 
short-lived 6^3/2 state of 133 Cs [f0(. Except for the case 
of the monovalent Cs atom, deducing octupole moments 
from measured HFS constants presents a formidable chal- 
lenge, owing to the fact that the prerequisite atomic- 
structure calculations become inaccurate for multi-valent 
atoms. As a result, previous analyses focused primarily 
on ratios of octupole moments for various isotopes of the 
same atom, because the atomic-structure factor cancels 
out when ratios of HFS constants are formed. By con- 
trast, the divalent alkaline-earth atoms considered here 
potentially yield direct values of the nuclear octupole mo- 
ments. 

3 i-2 states in alkaline- earth atoms are well suited for 
extracting nuclear magnetic octupole moments. We list 
the relevant isotopes in Table [IJ The atomic lifetimes are 
longer than 100 seconds and are long enough to allow for 
a sufficiently high spectral resolution. Moreover, success- 
ful trapping of the 3 P2 alkaline-earth atoms Sr and Ca 
has been recently reported in a number of laboratories 
and there are ongoing efforts with Mg [H M EE M- 
A typical trapping time is longer than 10 seconds. Since 
the magnetic trap perturbs the atomic level structure, 
a plausible experiment could involve microwave spec- 
troscopy of freely-falling atoms when trapping fields are 
turned off. Even then, the interrogation time is suffi- 
ciently long to permit a determination of the magnetic 
octupole HFS constant. To extract the nuclear octupole 
moment from the measured HFS constant, one requires 
electronic structure calculations. As shown in this paper, 
such supporting calculations can be carried out with an 
accuracy of a few per cent. 
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TABLE I: Nuclear parameters and lifetimes of the 3 P2 states of the stable isotopes of Be, Mg, Ca, and Sr. / is the nuclear 
spin; the superscript ir represents the parity. Experimental dipole and quadrupole values are from Ref. [H]. Octupole and 
hexadecapole values are from a single-particle approximation (see Appendix[S}. Lifetimes for the 3 p2 states are from Ref. 



Isotope 
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M cxp ' (Miv) 


Q oxp ' (b) 


n s p - (b x fi N ) 


n s p (b 2 ) 


3 P 2 Lifetime (s) 


9 Bc 


3/2" 


-1.177492(17) 


+0.053(3) 


-0.073 


0.00 




25 Mg 


5/2+ 


-0.85545(8) 


+0.201(3) 


-0.15 


0.00 


1.06 x 10 3 


43 Ca 


7/2- 


-1.317643(7) 


-0.049(5) 


-0.23 


0.00 


5.13 x 10 2 


87 Sr 


9/2+ 


-1.0936030(13) 


+0.335(20) 


-0.38 


0.00 


1.29 x 10 2 



This paper is organized as follows. In Section [IT] we 
recapitulate the major results from the relativistic the- 
ory of the hyperfine interaction. The consideration in- 
cludes first-order effects of nuclear dipole, quadrupole, 
octupole, and hexadecapole moments and second-order 
dipole-dipole and dipole-quadrupole effects. In Sec- 
tion IIIII we focus on the electronic structure aspects of 
the problem and evaluate electronic coupling constants 
using a configuration-interaction (CI) method coupled 
with many-body perturbation theory (MBPT). We re- 
fer to this approximation as the CI+MBPT method. 
Further, in Section IIV1 we discuss our results for the 
HFS of the 3 P 2 level in 87 Sr. Results for isotopes of Be, 
Mg, and Ca are tabulated in Appendix [U] Finally, con- 
clusions are drawn in Section |Vl The paper also has 
several appendices where we tabulate results of a tech- 
nical nature. Unless specified otherwise, atomic units, 
h = \e\ = m e = 1, and Gaussian electromagnetic units 
are employed throughout. 

II. NUCLEAR MOMENTS AND THE 
HYPERFINE INTERACTION 

The electric and vector potential of a nucleus modeled 
as a collection of point-like multipole moments may be 
expressed as 

Mr) = -iY,(- i r^ T] [^cZ(r)M k ^, (1) 

where Cfc./j and C k are normalized spherical harmonics 
and vector spherical harmonics, respectively (see, for ex- 
ample, Ref. [T3| ) - In these expressions, Q k ^ and M ktll 
are components of the nuclear electric and magnetic 2 fc - 
pole (MJ and EJ) moment operators, respectively. Each 
of the moments is an irreducible tensor operator of rank 
k. Parity and time-reversal symmetries require that k is 
even for the electric moments and k is odd for the mag- 
netic moments. Components of these tensors are conven- 
tionally parameterized using c-numbers that are defined 
as expectation values of the zero components of the op- 
erators in a nuclear stretched state \I, Mj = I). We will 
employ the following notation for the "stretched" matrix 



element of a tensor operator O k ^\ 

(O k )! = (I, Mi = I\O kfi \I, Mi = I). 

In particular, the magnetic-dipole, electric-quadrupole, 
magnetic-octupole, and electric-hexadecapole moments 
of the nucleus are defined as 

H = (Ali)/, 
Q = 2{Q 2 )i, 

n = -(M 3 )i, 
n = {Q 4 )i. 

The stretched matrix elements are related to the reduced 
matrix elements by 

(i\\o k \\i)= (-ioi) <°*>'- 

The interaction Hamiltonian for a single electron in an 
electromagnetic field is given by 

h em (r) = a ■ A (r) - ip (r) , 

and the total electromagnetic interaction Hamiltonian is 
then given by 

Hem ^ ^cm (^i) ) 

i 

where the sum runs over all electron coordinates. Sub- 
stituting the multipolar expansions |T|) into the electro- 
magnetic interaction Hamiltonian, we arrive at an ex- 
pression for rotationally-invariant hyperfine-interaction 
(HFI) Hamiltonian in the form of 

Here the spherical tensors (of rank k) T| act on elec- 
tronic coordinates and T£ operate in the nuclear space. 
We identify = M^'ior odd k and T™^ = e*,„ for 
even k. Explicitly, electronic tensors read 

^fc.ju — ^ tt,n ( r i) > 
i 

with 

!--pt+rC k ^ (f) electric (even fc), 

-pUy ^ a ■ C Z ( f ) ma § netic (° dd fc )- 

(2) 
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Matrix elements of these operators for Dirac bi-spinors 
are listed in Appendix [B] 

Now we recapitulate the application of perturbation 
theory to determining the hyperfine structure of atomic 
levels. In the presence of the multipolar fields produced 
by the nucleus, the total electronic angular momentum J 
is no longer conserved. The conserved angular momen- 
tum includes the nuclear angular momentum I; explicitly 
this conserved angular momentum is F = I + J. The 
proper basis consists of states \^IJFM F ) generated by 
coupling nuclear and electronic wave functions, 



\jIJFMf) = 



E 

Mj,Mi 



c 



FM F 

JMj-IM; 



\jJMj) \IMi) 



with 7 encapsulating remaining electronic quantum num- 
bers and the coupling coefficients being the conventional 
Clebsch- Gordon coefficients. For fixed values of J and /, 
the values of F range in | J — I\ < F < J + I, implying 
that an unperturbed level with angular momentum J is 
split into 2 J + 1 levels for J < I and into 21+1 levels 
otherwise. 

A matrix element of the HFI in the coupled basis is 



(yiJ'F'M' F \H HFI \jIJFM F ) 



->f'fom'm f 



x(-l) 



I+J+F 



E 



F J I 
k I J 1 



(yj , ||T fc e || 7 J)(7||T fe "||7), 
(3) 



the ^-symbols reflecting the scalar character of the HFI 
in the coupled basis. Hyperfine corrections to an un- 
perturbed level are given to first order by the diagonal 
matrix elements of -/7hfi- For convenience, we write this 
correction in terms of the i^-independent product of the 
stretched matrix elements 



W. 



(i) 



(-ylJFMp | //hpi | "fUFMp ) 



with 



X k {IJF) = (-1) 



I+J+F 





j 






i 


'A 



J k J 
~J J 



The first-order independent effects are conventionally 
parameterized in terms of the hyperfine structure con- 
stants A,B,C,D, - ■ ■ . Successive labels correspond to 
contributions of a multipolc of increasing multipolarity, 
e.g., A is due to the magnetic dipole moment, B due to 
the electric quadrupole moment, etc. The constants, up 
to D, are defined as in |l8[: 



A 
B 
C 
D 



l 

/./ ' 
4' 



(T?)i{T§)j 



J - 77 



1 

1.1 

2 



Q(T 2 e )j 

n(T|), 7 



(4) 



The second-order (in the HFI) correction for the state 
described by electronic quantum numbers 7 and J is 



W, 



,2) = I (j'IJ'FM F \Hj JFI \<yIJFM F ) 



E 1 j — E 1 >j> 



where the sum excludes the case (7' J') = ("fJ)- With 
Eq. this can be expressed explicitly in terms of re- 
duced matrix elements as 



(2) 



fci ,fc 2 

.,' T'l 





J 






J 




\ h 


I 


'A 


\ k 2 


I 


'A 



x(7 , J , ||^ 1 ll7J)<7'J'll^ 2 ||7J)</||^ 1 l|/)(/||^ 2 ||/). 



For the case of interest, where J = 2, the sum over (7' J') 
is dominated by contributions from the adjacent fine 
structure levels (7 1) and the sum over kik 2 is dominated 
by the dipole-dipole (fci = k 2 = 1) term. If we limit our- 
selves to the dipole-dipole and dipole-quadrupole terms, 
then we may express the second-order correction as 



W, 



(2) 



F J I 
1 I J -1 



V 



F J I 
1 I J- 1 



F J I 
2 I J — I 



where 77 and £ are -F-independent terms given by 

(J+1)(2J+1) 2 |<7J-1||TT|| 7 J}| 2 
V = 7 M 



E. 



E. 



7J-1 



c _ (J+l)(27+l) t / 2J + 3 „ 



fiQ 



I V 21 - 1 

(7J-i||rf|| 7 J)(7J-i||r|||7J) 

E 1 j — E 1 j-\ 



(5) 



(G) 



In the above expressions, 7 denotes a given fine structure 
term such as 3 P 2 and that J' = J — 1 = 1. 



III. CI+MBPT ELECTRONIC WAVE 
FUNCTIONS 

A precise analysis of the hyperfine structure depends 
on accurate electronic wave functions as well as nu- 
clear multipole moments. To obtain accurate electronic 
wave functions, we used a combined method of config- 
uration interaction (CI) and many-body perturbation 
theory (MBPT), which we refer to as CI+MBPT. The 
CI+MBPT method is described in detail, for example, 
111 Refs. [USQi. 

Here, we restrict the presentation to 
a qualitative discussion and then apply the CI+MBPT 
method to determination of the HFS couplings. 

The accuracy of the CI method is limited only by the 
completeness of the set of configurations used. In the 
context of CI+MBPT we refer to the sub-space contain- 
ing the configurations to be treated by the CI method 
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as the model space. In CI+MBPT, additional contri- 
butions from configurations which are not in the model 
space can be accounted for by MBPT. For our purposes, 
we are interested in finding accurate wave functions. It 
is worth noting that the wave functions determined by 
the CI+MBPT analysis lie completely within the model 
space. 

In deciding which configurations are to be included in 
the model space, we consider two things. First, configu- 
rations which interact strongly with the configuration of 
interest must be included in the model space. The main 
purpose here is to expand the model space to a degree 
in which the wave function can be described accurately. 
Second, configurations with energies which are nearly de- 
generate with the configuration of interest should also be 
included in the model space. These configurations lead 
to convergence problems in MBPT and are treated non- 
perturbatively by the CI method. Of course, "strongly" 
and "nearly degenerate" are relative terms which depend 
on the level of accuracy desired. 

For the alkaline earth systems considered in this paper, 
we start with a lowest-order description of our system as 
two valence electrons outside a closed core in a central 
field. Configurations involving different excitations of the 
valence electrons outside the closed core tend to interact 
strongly and have relatively close energies. These config- 
urations compose our model space. Additional contribu- 
tions from configurations involving excitations from core 
electrons are then accounted for within the framework of 
MBPT. 




(a) (b) 




(0 

FIG. 1: Sample Brueckner-Goldstone diagrams included in 
the effective CI+MBPT Hamiltonian. (a) Coulomb interac- 
tion between the two valence electrons. The double-arrows 
indicate that the single-particle states originate/terminate in 
the model space, (b) Brueckner (core-polarization) diagrams, 
(c) Screening diagrams. 

The relevant diagrams are shown in Fig. [TJ The 
Coulomb interaction between the two valence electrons is 
taken into account essentially to all orders of perturba- 
tion theory and the remaining many-body effects (involv- 
ing core excitations in the intermediate state) are treated 
in the second order of MBPT. In Fig.Q]the single-particle 
states are assumed to be generated in the frozen-core 



TABLE II: Energies of Sr I obtained from CI+MBPT are 
compared to experimental values. The table is partitioned 
into states of definite J and parity. Energies are referenced 
from the ground 5s 2 ^So state. All energy values are in cm -1 . 



J Parity 


Level 


CI+MBPT 


Expt. 


Diff. 


even 


5s6s 


L S 


30766 


30592 


174 




5p 2 


3 Po 


35798 


35193 


605 




5p 2 


'So 


37553 


37160 


393 


1 odd 


5s5p 


3 Pi 


14841 


14504 


337 




5s5p 


'Pi 


21818 


21699 


119 




5s6p 


s Pi 


34062 


33868 


194 




5s6p 


'Pi 


34293 


34098 


195 


2 odd 


5s5p 


3 P2 


15212 


14899 


313 




4d5p 




33836 


33267 


569 




4d5p 


'D 2 


34149 


33827 


322 



(V N ~ 2 ) Dirac-Hartree-Fock (DHF) approximation where 
the so-called one-body (loop) diagrams vanish identically. 

The specific approach that we used here is somewhat 
more sophisticated 20]: as the single-particle orbitals 
we employ the so-called Brueckner orbitals. This ap- 
proach incorporates the diagrams Fig. [ljb) from the on- 
set and also subsumes higher-order chains. Then the 
single-particle states (vertical lines) in the remaining di- 
agrams are replaced by the Brueckner orbitals. The 
entire scheme incorporates all the second-order MBPT 
effects and includes certain classes of diagrams to all 
orders. Technically, generating Brueckner orbitals re- 
quires an energy-dependent self-energy operator (see, 
e.g., Ref. [20]), S(e K ), which generally leads to non- 
orthogonal basis sets. We avoid this problem by fixing 
e K to the V N ~ 2 DHF energy of the lowest-energy valence 
orbital for a given angular symmetry n (k is defined ex- 
plicitly in Appendix [Bj) . 

To illustrate the predictive capabilities of the 
CI+MBPT method, in Table |TT] we compare the theo- 
retical energy values with the experimental values for Sr. 
The difference does not exceed a few 100 cm -1 for all 
the calculated energy levels. While the accuracy can be 
improved further, this level is sufficient for the goals of 
the present paper. 

Solving the CI+MBPT problem reduces to diagonaliz- 
ing the effective Hamiltonian in the model space. With 
the CI+MBPT wave functions we proceed to evaluat- 
ing matrix elements. Details are given in Ref. [2(j; the 
method builds upon the formalism originally developed 
for He- like systems [21]. We also include an important 
chain of diagrams of the random-phase approximation 
(RPA) in the evaluation of matrix elements. Qualita- 
tively, RPA accounts for screening of externally-applied 
fields (in our particular case these are the multipolar nu- 
clear fields) by the core electrons. 

Numerically, the calculations were carried out using B- 
spline basis sets. It is worth mentioning our modification 
to the original scheme [22j of generating the orbital sets. 
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In that scheme, boundary conditions are imposed on the 
small and large components of the wave function which 
were found to effectively dispense of spurious states in the 
orbital sets (more accurately, these spurious states were 
shifted towards the end of the spectrum). The disadvan- 
tage that follows is that the resulting wave functions are 
not highly accurate near the nucleus, leading to inaccu- 
racies in the HFI integrals of Appendix [B] To remedy 
this problem, we followed the prescription for creating a 
dual kinetic-balance (DKB) basis set, as introduced by 
Shabaev et. al. [23[. The DKB basis set is devoid of 
the problem of spurious states and capable of accurately 
representing the wave functions near the nucleus. The 
specific formulas used here for generating DHF orbital 
sets with a DKB basis set derived from P-splines are 
presented in Ref. [24|. 

Another technical issue is the appearance of "intruder" 
states in screening diagrams, Fig. [3Jc). Intermediate 
states (the diagram is cut across horizontally) there in- 
volve core-excited states. Since our model (CI) space is 
essentially complete, such core-excited states become em- 
bedded in the two-particle energy spectrum of the lowest- 
order model Hamiltonian. This leads to singular energy 
denominators in the diagrams Fig. |TJc). We remedied 
this problem by evaluating screening corrections only to 
the two particle states that had energies below the lowest- 
energy core excitation. For a typical basis set the result- 
ing subspace involved roughly 10% of the entire model 
space. Arguably, this prescription recovers most of the 
relevant correction since this low-energy subspace con- 
tains the dominant configurations. 



IV. HYPERFINE STRUCTURE OF B7 SR 

Using the techniques described in the previous sec- 
tions, we consider as a specific example the hyperfinc 
structure of the lowest-energy 3 P 2 level of 87 Sr. 



tion for 87 Sr in terms of these constants as 
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w- 
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65 r 






"21 C " 


—D 






3 



D, 



w< 



9/2 



+ Hl 

900 300 y/6 
32?7 4C \l 
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11/2 2 24 3 9 

13?7 



13C 



1100 550^6' 
W 13/2 = 9A+ ^B + C + D. 

The resulting HFS intervals 5Wp = Wp—Wp+i are given 
in terms of the hyperfine constants as 



6W< 



5/2 



-A H B C H D 

2 16 2 18 

in K 



900 300V6' 



5W 7/2 = - 9 -A+^B + ^C-^D 



14 2 

17/7 s l-i 



3300 220 V 2 



9/2 



SWu/2 



1 1 , 11 385 „ 

— A-\ CH D 

2 2 18 

V 



C 



13 

T 



24 



900 30V6' 
13 65 

3" C_ T 

13r7 13C 



-D 



1100 550V6' 



(7) 



Similar expressions for the hyperfine intervals in 3 P2 
states of stable isotopes of other alkaline-earth atoms are 
given in Appendix [Cl Solving Eqs.([7]) for the hyperfine 
constants C and D, one finds 



Extracting the magnetic-octupole constant from 
measurement of hyperfine intervals 



The nuclear spin of the stable 87 Sr isotope is I = 9/2, 
and so there are five hyperfine structure levels F = 
5/2, • • • , 13/2. First-order corrections are then charac- 
terized by four hyperfine structure constants A, B, C, and 
D. In addition, the second-order dipole-dipole interac- 
tion, characterized by r\ and the dipole-quadrupole inter- 
action, characterized by C> mix the 3 P 2 state with the 
nearby fine structure 3 Pi state. Using the expressions 
given m Section HH we may write the hyperfine correc- 



C 



D 



50 5/ 
21 



275 
3 

35V 



-8W< 



7 

550 
147 



7/2 



9/2 



6W< 



3575 
9 



SW n 



/2 



5/2 



550 



5W- 



7/2 



275 



-SW ( 



9/2 



3575 



1375^6 



5W 11/2 



(8) 



(9) 



Expressions for the HFS constants C and D in terms 
of hyperfine intervals for 3 P 2 states in other isotopes of 
the alkaline-earth atoms are given in Appendix [C] The 
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constant C depends only on the second-order dipole- 
quadrupole interference term £ while D is independent 
of both 77 and £. This proposition is independent of nu- 
clear spin /, as shown in Appendix [Dl 

To reiterate, measuring hyperfine intervals of the 3 P2 
level should allow one to deduce the magnetic-octupole 
HFS constant C, limited only by the knowledge of £. 
With the aid of calculations of the electronic structure 
factor presented below one may extract the nuclear mag- 
netic octupole moment of interest. 



B. Electronic structure factors 

From Eqs. (01 and ([5]) we see that the HFS constants 
may be written in terms of products of electronic ma- 
trix elements and nuclear multipole moments. We eval- 
uate the electronic matrix elements using the relativistic 
many-body method described in Section [TTTT We present 
the results of our calculations at various levels of approx- 
imation in Table IIIII 

In generating the DHF orbitals we included nuclear- 
size effects by assuming a Fermi charge distribution inside 
the nucleus. For computation of HFI integrals, however, 
we assumed a point-size nucleus. The observed effect on 
the choice of a particular model of nuclear distribution 
is below our theoretical error for solving the electronic- 
structure problem. 

There are several observations to be made with re- 
spect to the many-body calculations. First of all, we 
find that for Sr, the configuration 5p3/25si/2 provides 
the dominant (92%) contribution to the CI wave func- 
tion of the 3 P2 state. In particular, due to the angu- 
lar selection rules, the constants B and C are deter- 
mined by the matrix elements involving the 5j? 3 / 2 or- 
bital from this dominant configuration. By contrast, the 
dominant 5p 3 / 2 5si/2 configuration does not contribute to 
the electric-hexadecapole (E4) constant due to selection 
rules for single-particle matrix elements (see Appendix). 
Therefore the E4 electronic factor is strongly suppressed, 
as its value is accumulated entirely due to admixed con- 
figurations. 

We find that the BO (core-polarization) corrections 
universally increase the absolute value of all the con- 
stants. Qualitatively, core polarization describes an at- 
traction of the valence electron by the core. This attrac- 
tion leads to enhanced density closer to the nucleus and 
simultaneously larger hyperfine constants. Similarly, the 
computed RPA corrections show that the internal nu- 
clear fields are enhanced by virtual core excitations. The 
screening diagrams, Fig. QJc), qualitatively represent an 
effect of "cross-talking" between electrons via core polar- 
ization: a valence electron polarizes the core and this in- 
duced polarization attracts or repels another valence elec- 
tron. We see for the magnetic-dipole HFS constant that 
this effect is relatively weak compared to the other many- 
body corrections; however, its effect is more substantial 
for the electric-quadrupole and magnetic-octupole HFS 



constants. In all three of these cases, the screening con- 
tribution has the effect of decreasing the absolute value 
of the HFS constants. 

Finally, in Table IIIII we compare our ab initio results 
with experimental values for A and B. We find an 8% 
agreement for both constants. We believe that these ac- 
curacies are indicative of the theoretical error for the elec- 
tronic factor entering the magnetic-octupole constant C . 
The accuracy of computing the electronic factor for the 
HFS constant D is expected to be worse because the en- 
tire value is accumulated due to correlation effects. 

We have carried out similar many-body calculations 
for the parameters rj and C entering the second-order 
correction to the hyperfine constants. For 87 Sr, we find 
7] = 6.65 MHz and ( = 0.529 MHz. These second-order 
corrections are scaled to the experimental A and B coef- 
ficients and are accurate to 2%. The second-order dipolc- 
quadrupole contribution to C is 



AC( 87 Sr 



1375^/6 



= 1.10(2) kHz. 



(10) 



In this section we have described our calculations for 
87 Sr. The corresponding results for 9 Be, 25 Mg, and 43 Ca 
are given in Appendix [Cj 



V. DISCUSSION 

At this point we combine the computed electronic 
structure factors for the magnetic-octupole constants 
fTable UTIl and Table HV|) and the nuclear shell-model pre- 
diction for the M3 moment (Table [I] and Appendix [A"|). 
We find that 



C( 9 Be) = 


-3.57 x 10 


C( 25 Mg) = 


-2.57Hz, 


C( 43 Ca) = 


-16.2Hz, 


C( 87 Sr) = 


-201 Hz. 



(11) 



In particular for Sr, using Eq.([8]), we may deduce that an 
experimental sensitivity in measuring the HFS intervals 
on the order of a$w ~ 1 kHz would result in an uncer- 
tainty in the C constant on the order of o~c ~ 0.11 asw ~ 
100 Hz and would thus be capable of revealing the effects 
of a C constant of the predicted magnitude. 

The expression for the constant of interest C in 
term of HFS splittings contains the second-order dipole- 
quadrupole correction, see, e.g., Eq.© for Sr. If the 
experiment measures HFS intervals with a vanishingly 
small error bar, the extraction of would be limited by 
the error in this correction. Our estimated error bar for 
Sr, Eq. (fTTJ|) is 20 Hz, which translated into 10% error 
bar for C( 87 Sr) of predicted magnitude. Similar conclu- 
sion holds for Ca, while for Mg our estimated uncertainty 
in the dipole-quadrupole correction is comparable with 
the predicted size of C . For Be, our present uncertainty 
of 0.2 Hz in the dipole-quadrupole correction precludes 
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TABLE III: Breakdown of many-body corrections to hyperfme structure constants of 8r Sr 5s5p 3 P2 state. We used /i = 
— 1.0936/iiv and Q = 0.305(2) b (Ref. [25|]) in tabulating A and B constants. The final result is compared with experimental 
values from Ref. [2(|. CI-DHF corresponds to CI values computed using single-particle basis generated in the frozen-core 
(V N ~ 2 ) DHF potential. 





A, MHz 


B, MHz 


C/n, MHz/Vjv x b) 


D/IL, MHz/b 2 




CI-DHF 


-147.1 


35.6 


3.54 x 10" 4 


0.54 x 10" 


12 








Many-body corrections 






A BO 


-41.9 


10.2 


1.03 x 10" 4 


0.85 x 10~ 


12 


A Screen 


0.4 


-4.6 


-0.42 x 10~ 4 


2.71 x 10" 


12 


A RPA 


-42.0 


21.0 


1.15 x 10~ 4 


0.55 x 10 - 


12 


Final 


-230.6 


62.2 


5.30 x 10~ 4 


4.65 x 10" 


12 


Experiment 


-212.765(1) 


67.215(15) 









clean extraction of the octupole moment. Notice, how- 
ever, that the nuclear shell model estimates of the nu- 
clear octupole moment may be unreliable. If the 133 Cs 
experiment [lfj is of any indication, the "true" size of 
the octupole constant may be much larger (factor of 40) 
than predicted. Then the dipole-quadrupole corrections 
become mostly irrelevant. 

We emphasize that the values (fTTj) are only estimates 
based on the nuclear-shell model; measuring C would 
show deviations from these estimates. In a particular 
case of 133 Cs the measured and the predicted values dif- 
fered by a factor of 40 [l(| • It remains to be seen if such 
large deviations from the nuclear shell model would be 
revealed experimentally for the nuclei considered in the 
present work. 
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APPENDIX A: NUCLEAR MOMENTS FROM A 
SINGLE PARTICLE MODEL 



A crude approximation to the nuclear moments can be 
achieved by representing the nucleus by a single nucleon. 
For even-odd (even number of protons, odd number of 
neutrons) nuclei we use a single neutron, and for odd- 
even nuclei we use a single proton. Using formulas by 
Schwartz (27| . we may write the moments in this single 



particle model as 



= UnI x 



f'P- = -e(r 2 )gi 




90/(2/ + 2)] 



for 7 



for 1 = 1-1, 



2 (27 + 4) (27 + 2) 
f(/+2)[(I-§)fl,+5.] 

\(I-l)[(I + l) gi -g s ] 



for 7 
for 7 



I - i 

1 2 ' 



jjs.p. 



3 (27- l)(27-3) 
8 (27 + 4) (2/ + 2)' 

where gi = +1, g s = 5.58 for a proton and gi = 0, 
g s = —3.83 for a neutron. All of the stable isotopes 
considered in this paper have even-odd nuclei. This has 
the immediate consequence Q s p - = and IP p = for 
all isotopes. Furthermore, an examination of the momen- 
tum 7 and parity ir in Table [J reveals that the nucleon 
for each isotope must have an orbital momentum I sat- 
isfying 1 = 1 + 1/2. With this additional consideration, 
/i s p ' = — 1.92/ijv for all isotopes, and the expression for 
the octuple moment is reduced to 

Approximating the root-mean-square value of the nuclear 
radii, (r 2 ) 1 / 2 , as (in units of fm) 2.52, 3.05, 3.48, and 4.24 
for the cases of 9 Be, 25 Mg, 43 Ca, and 87 Sr, respectively, 
yields the values for J7 s p - given in Table fl] 

In the single particle model the electromagnetic mo- 
ments of the nuclei are given by the appropriate ex- 
pectation values for the valence nucleon shell. In this 
model the even-odd nuclei would have electromagnetic 
moments determined by the valence neutron. In particu- 
lar, since the neutron doesn't have an electric charge, all 
electric moments vanish. Certainly, the observed nonzero 
Q-moment provides an information on such quantities 
as nuclear deformation. Similarly, an observation of the 
electric hexadecapole moment, which is zero in the single 
particle approximation, should provide similar informa- 
tion on the nuclear distortion. 
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APPENDIX B: MATRIX ELEMENTS OF THE 
ELECTRONIC TENSOR OPERATOR 



2:. 



Mg, / = 5/2: 



Here we compile expressions for the matrix elements 
of the single-particle electronic HFI coupling operators 
t% (r) given in Eq. We use the conventional 

parametrization of the Dirac bi-spinors, 



mum) 



1 ( iP nK {r) Km (f) 
r I Qn K (r) fl- Km (f) 



where k — (I — j) (2j + 1) and Q are spinor functions. 
With this parametrization, we find that the reduced ma- 
trix elements for the HFI couplings to electric moments 
of the nucleus are given by 

(nV||^|K) = -W\\C k \\ K ) 

dr 

\Pn' k' Ptlk ~\~ Qn' n'Qnn) 3 

and for couplings to magnetic moments these are 



SW 1/2 = -^A 



SW- 



3/2 



sw 5/2 = 



SW 7/2 = 



- 7 -A- 
2 



9 54 

T) 



c 



-A + -B — 3C — 75D 
2 2 



100 50V2' 

13r? C 



7 49 



1260 210V2' 

V IK 
900 300 V2' 



27 



27 



— B C -9D 



40 



3r? 33C 
140 140\/2' 



<nV||tg|M = (/c'llCfcH-zs) 



Selection rules for these matrix elements follow from 
those for the matrix elements of the C-tensor: \j — j'\ < 
k < j + j' and the sum I + V must be even. 



c = - To sw l/2 



D = -SW 1/2 



70 



-SW: 



3/2 



20 

5 

252 v 



-SWs 



5/2 



SW* 



490 



-SW; 



3/2 



980 
1 

1764' 



7/2 



SW5/2 



-SW: 



7/2 



350V2 



APPENDIX C: HYPERFINE STRUCTURE OF 
THE 3 P 2 STATE FOR 9 BE, 25 MG, AND 43 CA 

This appendix contains a compilation of expressions 
relating the hyperfine intervals and the hyperfine struc- 
ture constants for the 3 P 2 states of 9 Be, 25 Mg, and 43 Ca. 
Computed HFS constants for these isotopes, given in Ta- 
ble |lVl are also included. 

9 Bc, / = 3/2: 



SW 1/2 



SW 3/2 = 



-A- 



7 



B - 28C 



-B + 20C- 



200 ' 



600 
V 



SW, 



7 



5/2 



A 



7 



B-7C 



120 40 ' 
7y 7C 
200 200V3 1 



C = -—5W 1/2 + — SW 3/2 - — SW 5/2 + — 2_ 
50 1/2 50 3/2 175 5/2 500^3 

For 9 Be, we find that -q = 25.09(4) MHz and ( = 
0.2939(2) MHz, leading to a value AC = 0.3394(2) kHz 
for the second-order correction to C. 



For 25 Mg, we find that 77 = 5.37(1) MHz and C = 
0.333(1) MHz, leading to a value AC = 0.671(2) kHz 
for the second-order correction to C. 



43 



Ca, / = 7/2: 



SW- 



3/2 



SW, 



5/2 



SW- 



7/2 



-A - 



-A 



-A 



25 b 55 r 



275 



-D 



-2- + -^- - 

112 112V 3' 



-B - AAD 



c 



144 48\/l5' 



56 



-B 



7 7 



V 



13C. 3 



SW, 



33 



9/2 



'Ia-^b 

2 56 



1680 560 V 5' 

7 7 
H77 11c 



720 240VT5 



9 



117 

20 1 15 1 220 1 



C 



120V15 



660 



SW< 



9/2 



For 43 Ca, we find that r) = 8.43(3) MHz and C = 
0.085(1) MHz, leading to a value AC = -0.183(3) kHz 
for the second-order correction to C. 



The second term becomes 





J 






j 






I 


'A 




i 


A 



The third term becomes 



*?(-!) 



2{I+J)+k+l 



?k,k' 



£(-*) 



(2k + 1) (2k + 1) 



F-I-J-k-1 



(2F+1) 



APPENDIX D: PROOF THAT HFS CONSTANTS 
C AND D MAY BE DEFINED COMPLETELY IN 
TERMS OF THE HFS INTERVALS 

In this section we prove that the HFS constants C 
and D can be expressed uniquely in terms of the HFS 
intervals even when the second-order dipole-dipole fine 
structure term may not be neglected. This is a nontrivial 
statement, as the number of linear equations for HFS 
intervals is less than the number of fitting parameters. 
For example, for 87 Sr we find that there are four HFS 
intervals expressed in terms of five fitting parameters (see 
Eq.©). 

We start by defining our HFS levels from a new en- 
ergy offset, W'p = Wf + A. The constant A is arbitrary 
and its particular choice will be shown not to affect the 
conclusions; consequently, a knowledge of the HFS level 
intervals with a convenient choice of A is sufficient to 
completely define W'p for all F. Including all first-order 
terms as well as the second-order dipole-dipole fine struc- 
ture term, the levels W' F can be written as 



WL 



A + (-1 



J+J+F 



E 



F J I 
k' I J 



F J I 
1 I J —I 



V, 



(Dl) 



where Z k = ^J\\T^\\ 7 J)(I\\T^\\I). From Eq. © we 
see that Z\ is proportional to A, Z^ is proportional to 
B, etc.. The next step is to multiply every term in Eq. 

- by (-1) /+J+F (2F + l){l J T \ ), with k f 0, 



v k 1 J , 

and sum over all -F-values. Here we analyze the effect 
of this procedure on the individual terms of the right 
hand side of Eq. (|D1[) ; to do so, we incorporate various 
well-known sum rules of six-j symbols. The first term 
becomes 



A-EM 



J+J+F 



(2F + 1) 



F J I 
k I J 



A • 5 kt0 yJ(2I+l) (2J+1) = 0. 



F J 



F J I 
1 I J -I 



1 I J- 1 j 

^2(/+.7) + fc+l / 1 1 k 



r,(-iy 



J J J -I 



F J I 
k I J 

1 1 k 
III 



The resulting equation may then be solved for Z k , giving 



Z k = (2k + l)^(-l) I+J+F (2F+l)[ F k i j}w' F 



+ (-1) 



2(/+J) + fc 



(2fc+l) 



1 1 k 
J J J-l 



1 1 k 
III 



First, we note that this expression does not depend on the 
specific choice of A. Second, we note that for the case of 
k > 2, the triangular condition is not satisfied along the 
top rows of the last four six-j symbols in the last expres- 
sion. Since these six-j symbols are then equal to zero, Z k 
is completely defined by the values of W' F . Equivalently, 
we may conclude that the HFS constants C D. ... may be 
expressed completely in terms of the HFS intervals, and 
these are the same expressions that would be obtained 
regardless of the inclusion of rj. A more general proof 
can easily be given to show that with the inclusion of 
second-order dipole-quadrupole terms, C can no longer 
be expressed completely in terms of the intervals, while 
the expression for D in terms of the intervals would still 
remain valid (and so on to higher second-order terms if 
desired). 

The above conclusion can also be drawn from the for- 
mulation of second-order effects as in Ref . [3l| , in which 
the second-order effects are used to describe the differ- 
ence between a measured value of a HFS constant (based 
on first-order perturbation theory) compared to its actual 
value. This already assumes that all measurable HFI ef- 
fects are completely described by first- and second-order 
perturbation theory. Further assuming that only second- 
order terms of the dipole-dipole type contribute to mea- 
surable effects shows that measured HFS constants differ 
from actual HFS constants only for the cases of constants 
A and B and not for higher constants. 
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TABLE IV: Theoretical and experimental hyperfine structure constants for the 3 P2 states of 9 Be, 25 Mg, and 43 Ca. Theoretical 
values include the many-body effects discussed in Section ITTT1 







A, MHz 


B, MHz 


C/fi, MHz/ (/Mr x b) 


D/n, MHz/b 2 


9 Be 


Tlisory 
Expt. a 


— 119.7 

-124.5368(17) 


1.43 

1.429(8) 


4 89 x 10 _Y 




25 Mg 


Theory 
Expt. 6 


-127.5 
-128.445(5) 


15.8 

16.009(5) 


1.71 x 10~ 5 


-1.39 x 10~ 14 


43 Ca 


Theory 
Expt." 


-179.9 
-171.962(2) 


-5.50 
-5.436(8) 


7.03 x 10" 5 


7.83 x 10~ 13 



a Ref. Q 
6 Ref. 2EJ 
c Ref. [13 
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